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Scattering processes in antiprotonic hydrogen - hydrogen atom collisionl^ 


V.P. Popov and V.N. Pomerantsev 
Institute of Nuclear Physics, Moscow State University 

The elastic scattering, Stark transitions and Coulomb deexcitation of excited antiprotonic hy¬ 
drogen atom in collisions with hydrogenic atom have been studied in the framework of the fully 
quantum-mechanical close-coupling method for the first time. The total cross sections (E) 

and averaged on the initial angular momentum I cross sections (E) have been calculated for 

the initial states of {pp)n atoms with the principal quantum number n = 3 — 14 and at the relative 
energies E — 0.05 — 50 eV. The energy shifts of the ns states due to the strong interaction and 
relativistic effects are taken into account. Some of our results are compared with the semiclassical 
calculations. 


I. INTRODUCTION 

The slowing down and Coulomb capture of the negative particle M~ {M~ = p,~, 7r“, etc.) in hydrogen 
media lead to the formation of the M“-molecular complex the decay of which results in the exotic atom 
formation in highly excited states with the principal quantum number n ~ where p is the reduced 
mass of the exotic atom. Their initial nZ-population and kinetic energy distribution of the exotic atom 
are defined by the competition of different decay modes of this complex. The further destiny of the exotic 
atom depends on the kinetics of the processes occurring in the deexcitation cascade. The experimental 
data are mainly appropriate to the last stage of the atomic cascade, such as X-ray yields and the products 
of the weak or strong interaction of the exotic particle in the low angular momentum states with hydrogen 
isotopes. 

Hadronic hydrogen atoms are of special interest among exotic atoms because they have the simplest 
structure and are the probe in the investigations of the various aspects of both the exotic atom physics 
and the elementary hadron-nucleon interactions at zero energy. In particular, in order to analyze the 
precision spectroscopy experimental data 0 the kinetic energy distribution must be taken into account. 
The velocity of the exotic hydrogen atom plays an important role due to the effect of the Stark transitions 
on the L X-rays yields and the Doppler broadening of the L-lines owing to the preceding Coulomb 
deexcitation transitions. The energy release in the last process leads to an acceleration of colliding 
partners. So the reliable theoretical backgrounds on the processes both in low-lying and in highly excited 
states are required for the detailed and proper analysis of these data. 

In this paper we present the first step to ab initio quantum-mechanical treatment of non-reactive 
scattering processes of the excited antiprotonic hydrogen atom in collisions with the hydrogenic atom in 
the ground state: 

- elastic scattering 


{ax)ni + {be )is ^ {ax)ni + {be )is; (1) 

- Stark mixing 

{ax)ni + {be~)is {ax)nv + {be~)is', (2) 

- Coulomb deexcitation 

{ax)ni + {be~)is {ax)n'i' + {be~)is. (3) 

Here (a, 6) = {p,d,t) are hydrogen isotopes and x = K~,p] {n,l) and Is are the principal and orbital 
quantum numbers of the exotic and hydrogenic atom, respectively. The processes (1) - (2) decelerate and 
accelerate while the Coulomb deexcitation (3) accelerates the exotic atoms, influencing their quantum 
number and energy distributions during the cascade. The last process has attracted particular attention 
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especially after the ’’hot” irp atoms with the kinetic energy up to 200 eV were found experimentally Q- 
Due to the similarity of the general features of the exotic atoms the Coulomb deexcitation process must 
be also taken into account for the other exotic atoms. 

Starting from the classical paper by Leon and Bethe [3], Stark transitions have been treated in the semi- 
classical straight-line trajectory approximation (see [4] and references therein). The first fully quantum- 
mechanical treatment of the processes (1) - (2) based on the adiabatic description was given in [5-8]. 
Recently!^ ^ the elastic scattering and Stark transitions (for n = 2 — 5) have also been studied in a 
close-coupling approach treating the interaction of the exotic hydrogen atom with the hydrogenic one in 
the dipole approximation with electron screenin g ta ken into account by the model. As for higher exotic 
atom states (n > 5), the semiclassical approach |l3 is used for the description of these processes. 

As concerning the acceleration process (3) in the muonic and hadronic hydrogen atoms, the parame¬ 
terization based on the calculations in the semiclassical model (see also [13]) is used for the low-lying 
states (n = 3 — 7) and the results of the classical-trajectory Monte Carlo model mill are used for 
higher exotic atom states in the cascade calculations. 

The main aim of this paper is to obtain the cross sections of the processes (l)-(3) for the excited 
antiprotonic atom beginning from the low collision energies in the framework of the fully quantum- 
mechanical approach. For this purpose the unified treatment of the elastic scattering, Stark transitions 
and Coulomb deexcitation within the close-coupling method has been used. This approach has been 
recently applied for the study of the differential and total cross sections of the elastic scattering. Stark 
transitions and Coulomb deexcitation in the collisions of the excited muonic 0 and pionic 0 hydrogen 
atoms with the hydrogen ones. In the following Section we briefly describe the close-coupling formalism. 
The results of the close-coupling calculations concerning the total cross sections of the processes (l)-(3) 
are presented and discussed in Section III. Finally, summary and concluding remarks are given in Section 
IV. 


II. CLOSE-COUPLING APPROACH 


A. Total wave function of the binary system in terms of basis states 


The total wave function 'I'(p, r, R) of the four-body system {ap + be~) satisfies the time independent 
Schrodinger equation with the Hamiltonian which after separating the center of mass motion can be 
written as 

H =-■;^Ar- f/ip(p)-f/le(i’) + R(i’, P, R-) (4) 

2m 

(m is the reduced mass of the system). Here we use the set of Jacobi coordinates (R, p, r): 

R Rr Rpa; P ^a: ^ ^b: 

where r^, r{,, r^, Tg are the radius-vectors of the nuclei, muon and electron in the lab-system and Rr, Rpa 
are the center of mass radius-vectors of hydrogen and exotic atoms, respectively. The Hamiltonians hp 
and he of the free exotic and hydrogen atoms, respectively, satisfy the Schrodinger equations 

hp*^nlm{p') — ^nl*^nlm{p'): (5) 

hePnjA^) = ene‘Pneie(r), (6) 

where ^nim{p) and y)„^(r) are the hydrogen-like wave functions of the exotic atom and hydrogen atom 
bound states, £ni and are the corresponding eigenvalues. In the present study £ni includes beyond 
the standard non-relativistic two-body Coulomb problem the energy shifts due to the strong interaction, 
vacuum polarization and finite size. It is worthwhile noting that in order to treat the hadronic ns states 
as normal asymptotic states in the scattering problem we take into consideration only the real part of 
the complex strong interaction energy shift. 

The interaction potential 


V(r, p, R) = Vab + Vpb + Rae + Vpe (7) 

includes the two-body Coulomb interactions between the particles from two colliding subsystems: 

Vab = — = \Ii+l'p-l^er\~^,Vpb = - — = -|R- Cp- (8) 

Tab rpb 

Vpe = — = |R - ^p-f $el'r\Ke = - — = - |R -f l^p J-^gr |\ 

ae 


( 9 ) 
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where the following notations are used: 

V = mplimp + nia), ^ = nia/(mp + rUa), ( 10 ) 

= rrielime +mb), = rribfime + rub), (11) 

{ma,mb,rnp and rUe are the masses of hydrogen isotopes, antiproton and electron, respectively). Atomic 
units (a.u) h = e = memb/{me + mb) = 1 will be used throughout the paper unless otherwise stated. 

In this paper, as well as in the previous studies [11, 12-15], we assume that the state of the target electron 
is fixed during the collision. The electron excitations can be taken into account in a straightforward 
manner. In a space-fixed coordinate frame we built the basis states from the eigenvectors of the operators 
he,hp,l‘^,L‘^,3‘^,Jz and the total parity tt with the corresponding eigenvalues Sis,Sn,l{l + l),L{L -|- 
1), J(J -|- l),M and (—1)^+^, respectively: 


\ls,nl,L: JM) = -^RUr)Rni{p)ytj^{p.tl), (12) 

V47r 

where 

J^//^(p,R) = Y,{lmL\\JM)Yi^{p)YLx{tC). (13) 

mX 

Here the orbital angular momentum 1 of {ap)ni is coupled with the orbital momentum L of the relative 
motion to give the total angular momentum, J = 1 -|- L. The explicit form of the radial hydrogen-like 
wave functions Rniif) will be given below. 

Then, for the fixed values of JjAIjir = (—1)^+^ the exact solution of the Schrodinger equation 

{E-H)^i^^r,p,K)=0, (14) 


is expanded as follows 


vI'i^-(r,p,R) = i^G;(r^(i?)|ls,nZ,L: JM), (15) 

nlL 

where the G'^ip{R) are the radial functions of the relative motion and the sum is restricted by the [I, L) 
values to satisfy the total parity conservation. This expansion leads to the coupled radial scattering 
equations 

GifLiR) =2mJ2 W;^'i'L'MLiR) Gi^L,{R), (16) 

^ ' n'l'L' 

where = 2m{Ecm + £niii — Sni) specifies the channel wave numbers; Ecm and are relative motion 
energy and exotic atom bound energy in the entrance channel, respectively. 

The radial functions (i?) satisfy the usual plane-wave boundary conditions at i? ^ 0 

Gi>L<(0)=0(~R^+i) (17) 


and at asymptotic distances {R oo) 

1 


G'^E,n'V L'{R) 


\Af 


- S\nl,L^ n'l', 


(18) 


where ki, kf are the wave numbers of initial and final channels and S'’[nl, L n'V, L') is the scattering 
matrix in the total angular momentum representation. Here and below the indices of the entrance channel 
and target electron state are omitted for brevity. 
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B. Potential matrix 


Here we present the derivation of the exact matrix of the interaction potentials involved in the close¬ 
coupling calculations. The interaction potential matrix coupling the asymptotic initial {nlL; J) 

and hnal {n'VL'\ J) channels is dehned by 


Wn'i'L,niL{R) = ^ J dr d /3 dRi?J,(p) 

X R) vir,p, R) ip, R), 

where the radial hydrogen-like wave functions are given explicitly by 

I n—l — 1 


Rniip) = exp(-p/na) ^ 


(a is the Bohr’ radius of the exotic atom in a.u.) with 

Nnl = 

and 

S,{n,l) = i-y 


/ 2 

^ — 1 

1 

1 

+ 

1 _ 

\na J 

2n 


1/2 


1 


(19) 


( 20 ) 


q\in - I - I- g)!(2/ -I- 1 -I- g)!' 

Averaging V (r, p, R) over Is-state of hydrogen atom leads to 

H(R,p) = —J^ dri?L(r)V^(r,p,R) = 

= - U-^^^^iR,p)} -{t/i.,j.JR,p) - t7_j,^^(R,p)}. 

se 


( 21 ) 


( 22 ) 


(23) 


Then we use the transformation 


p) 


(1 


f3 


|R -l- 0!p\ 


)e’ 


2|R.+ <»p| 


lim 

x^l 



2x|R.+ ctp| 


/3 


e P 
|R -l- oip\ 


(24) 


which allows us to apply the additional theorem for the spherical Bessel functions [it 


g—A|Ri-|-ri| 

|Rl + ri| 


47r 


^(-l)X(Ri)d^*r(fi)x 


X |Art+i/ 2 (Ai?i)/t+i/ 2 (Ari) +/t+i/ 2 (Ai?i) A:t+i/ 2 (Ari)|^^^^^| 


(25) 


{Ip{x) and Kpix) are the modihed spherical Bessel functions of the first and third kind). Furthermore, 
by substituting the Eqs.(20)-(25) into Eq.(19) we can integrate over the angular variables (R, p). Finally, 
applying the angular momentum algebra and integrating over p, we obtain: 


t=0 ^ 

X [(-l)*VFi(i?,Z2,^e;n^,R'^') - “ 

- - [(-l)*lFt(i?, Ve; nl, n'l') - WtiR, C, Ve] nl, n'l')] 1 

J 


( 26 ) 
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{tm is the maximum value of the allowed multipoles). Here the next notations are used: 


n-i-i / 2„' \rmn-i-i / 2n \ 

Wt{R,a,P;nl,n'l') = J\fni,n'i' 'Y] SmA'nJ) i —;— 7) V' Sm2in',l')[ —;— 7) 

V n + n' / V n + n' / 


/ \ mi n' —I' — 1 


n + n'J ^ ' \n + n' 

m\—0 ^ ^ m 2=0 ^ 

X {Hi(a:) A(n, n', a, /?)) — ht{x)jY (x, A(n, n', a, /?))+ 


+ Ft{x)Jl''\x, A(n, n', a, /?)) +/t(x) J 4 ’®(x, A(n, n', a, /?))} 
where a; = 2R/j3, s = I + V + mi + m 2 , L = 2L + 1; 




n + n 


(27) 


1 / O.^/ \ \ ^ "1"^ 


A(n, n', a, /?) = 


2nn' aa 
n + n' f3 ' 


(29) 


i7t(x) = (1 - 2t)ht(x) +xht+i{x)-, 


(30) 



Ft{x) = (1 - 2.t)ft{x) - xft+i{x). 

(31) 

The functions ht{x) 

and ft(x) are given by 



ht{x) = J —Kt+i/2{x) 

V 7 TX 

(32) 

and 

h{x) = \I^A+i/ 2 ix). 

(33) 


The radial integrals JY{x, A) are defined as follows: 


j‘’®(x,A)= [ yftiXy)dy, (34) 

j‘’*(x,A) =AJ4‘+'’*+'(a:,A), (35) 


pOO 

JY{x,X)= y''+'^e~yht{Xy)dy, 

J x/X 


(36) 


j*''*(x,A) = AJ3+^’"+^(a;,A) 

and calculated analytically using the power series for the modified Bessel functions. 


(37) 


C. Cross sections 


The transition amplitude from the initial state |n/m > to the final state > of the exotic atom 

can be defined by 

f{nlm^n'l'm'\]^i,'kf) = —^Y== 'Y A ~^{lmLX\JM){l'm'L'X'\JM)x 

JMWXX' 

X YlA^,)YL-y{^f)TAnlL ^ n'l'L'). 


( 38 ) 
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Here, ki and kf are the center of mass relative momenta in the initial and final channels; ki and kj are their 
unit vectors in the space-fixed system, respectively, and, finally, the transition matrix T'^(nlL —> n'VL') 
used here is given by 

L n'l', L') = Snn'SwSLL'Smm'Sx\i - S'^{nl, L n'V, L'). (39) 

In terms of the scattering amplitude (38) defined above, all the types of both the differential and total 
cross sections of the all processes under consideration for the transition from the initial {nl) state to the 
final {n'V) state are defined as: 
differential cross sections 

mm' 

partial cross sections 

^ni^n'AE) = (41) 

and the total cross sections for the nl n'V transition are obtained by summing the corresponding 
partial cross sections over the total angular momentum J: 


0'nl^n'l'{E) — ^ ^ CTnl^n'l' (E) ■ 
J 


(42) 


Finally, the averaged over the initial orbital angular momentum I cross sections for the n ^ n' transi¬ 
tions are then defined by summing over V and I {I = V for the elastic scattering and I ^ V for the Stark 
transitions) with the statistical weight (21 + l)/n^ in the case of the degenerated exotic atom states and 
with the weight (21 -|- l)/(n^ — 1) in the case when the energy shift of the ns state is taken into account: 


and 




Y, {2J+l)\T^{nlL^n'VL')\^, 

l' LL'J 


(43) 


J>0 


{E) = 


1 


kl 


^2 _ 1 

/> 0 ,1 ’ LL'J 


{2J+l)\T\nlL^n'VL')\^ 


(44) 


respectively. 


III. RESULTS 

The close-coupling method described in the previous Section has been used to obtain the total cross 
sections for the collisions of the pp atoms in excited states with hydrogen atoms. The present calculations 
had at least two goals: first, to apply the fully quantum-mechanical approach for the study of the processes 
(1) - (3) and, second, to clear the effect of the energy shifts of the ns states of the antiprotonic hydrogen 
atom on the cross sections of these processes. 

The coupled differential equations (16) are solved numerically by the Numerov method with the 
standing-wave boundary conditions involving the real it'-matrix. The corresponding T-matrix are ob¬ 
tained from itT-matrix using the matrix equation T = 2iK{I — iK)~^. In the calculations both the exact 
interaction matrix and all the open channels with n' < n have been taken into account. The closed 
channels are not considered in the present study. The close-coupling calculations have been carried out 
for the relative collision energies Ecm from 0.05 up to 50 eV and for the excited states with n = 3 -h 14. 
At all energies the convergence of the partial wave expansion was achieved and all the cross sections were 
calculated with the accuracy better than 0.1%. 

The results of the calculations are presented in Figures 1-4. In Fig. 1 we introduced the calculated total 
cross sections of the nl nV transitions for n = 8 at the kinetic energy = 1.4 eV both with and 
without the ns state energy shifts. The following measured world-average value for the spin-averaged 
shift 


ei^ = (721 ± 14)eV 
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FIG. 1: The total cross sections a^i^ni' for the collisions of the pp atom in the n = 8 state with hydrogen atom 
at iJcm = 1.4 eV. The dashed and dotted lines connect the points corresponding to the calculations both with and 
without taking into account the ns-state energy shifts, respectively. The dotted lines denote the results obtained 
without including the ns-states into the basis set. 



FIG. 2: The i-averaged cross sections of the elastic scattering for the collisions of the pp atom with the hydrogenic 
atom for the different values of the principal quantum number n 


was used in the present calculations and the energy shifts of the ns states are defined by txsjv?. 

Contrary to the (7r“p)-atom, the energy shifts of the ns states in (pp)-atom are repulsive, hence, the 
nl ns transitions are closed below the corresponding threshold and, besides, according to the present 
study (e.g., see Fig. 1) the Stark transitions both from the ns-states and to the ns-state at the same 
collisional energy are strongly suppressed at the kinetic energies above the threshold. The similar effect 
is also observed for the elastic np —>■ np transitions. The other transitions are practically unchangeable 
at fixed energy. The analog of this effect can be modeled by excluding ns-states from the basis states. 
At energy above a few thresholds the effect is much weaker. Therefore, the strong interaction effect in 















the antiprotonic hydrogen (the similar effect must be also observed in the kaonic hydrogen atom) results 
in the essential and important difference of the Stark transition role in the absorption from the ns states 
in contrast to pionic hydrogen atoms. The influence of the strong interaction shift enhances for the lower 
states and becomes less pronounced for the highly excited states of the antiprotonic atom. 



FIG. 3: The /-averaged cross sections of the Stark transitions for the collisions of the pp atom with the hydrogenic 
one. The results of the calculations m in the semiclassical model are shown for n = 8 with triangles 

In Figures 2 and 3 the energy dependence of the /-averaged total elastic and Stark cross sections 
are shown for the different principle quantum number values from n = 3upton = 12. Since the 
relative contribution of the ns state in the /-averaged cross sections are small, the calculations both with 
and without the energy shift are practically indistinguishable at the energies above the corresponding 
threshold. So, in Figs. 2 and 3 the small energy region (below the corresponding threshold) corresponds 
to the calculations without the energy shift taken into account. In Fig. 3 the /-averaged cross section for 
n=8 are compared with the calculations in the framework of the semiclassical model M- As a whole a 
fair agreement is observed, but the semiclassical model results in a different energy dependence especially 
at low energy collision and gives smaller cross sections than those obtained in the present approach. 

The dependence of the Coulomb deexcitation cross sections on energy obtained in the present study 
is illustrated in Fig. 4 for n = 8 and the different values of An =1, 2 and 3. The special features of 
these cross sections are the following: the similar energy dependence but sharper than that of the elastic 
scattering and Stark transitions (see also in Fig. 4); the contribution of the transitions with An > 1 
is comparable with the one for An =1 and approximately equal about 50%. The effect of the ns state 
energy shifts in the /-averaged Coulomb deexcitation cross sections are small for the same reason as it 
was discussed above (due to small statistical weight of the ns-state). In Fig. 4 we also compare our 
results with those obtained in the semiclassical model (we use the parameterization suggested in [13] 
which gives a fair description of the Coulomb cross sections from [11]) for the An =1 transition. The 
satisfactory agreement is observed, but this agreement is quite occasional and takes no place for other n 
values. The distribution over the final states n' is completely different from the semiclassical results [11] 
as it was illustrated in Fig. 4. The present calculations predict that An > 1 transitions give a substantial 
contribution to the Coulomb deexcitation of the highly excited antiprotonic hydrogen atom that is in 
agreement with our previous results for the muonic and pionic hydrogen [14,15] atoms. 


IV. CONCLUSION 


The unified treatment of the elastic scattering. Stark transitions and Coulomb deexcitation is presented 
in ab initio quantum-mechanical approach and for the first time the cross sections of these processes have 
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FIG. 4: The Z-averaged cross sections of Coulomb deexcitation with An = 1,2,3 for the collisions of the pp 
atom (n = 8) with the hydrogenic atom. The dashed line shows the £t used in cascade calculations fl^ for 
the transitions with An = land based on the mass-scaling of the results for the muonic atom. The present 
calculations of the Z-averaged elastic and Stark cross sections are shown for comparison 


been calculated for the highly excited antiprotonic hydrogen atom. The influence of the energy shifts of 
the ns states on these processes has been studied. We have found that strong interaction shifts in the 
antiprotonic hydrogen atom lead to substantial suppression of both the ns nV 0 and nV 0 —*■ ns 
transitions. At the same time the cross sections of the elastic scattering and Stark transitions for the the 
states with I > 2 are practically unchangeable. The present study is the first step to achieving a reliable 
theoretical input for realistic kinetics calculations. 

We are grateful to Prof. L.Ponomarev for the stimulating interest and support of this investigation. 
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